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ABSTRACT 

A  global,  six-level,  primitive  equatiDn  model  was  numer- 
ically  integrated  to  examine  the  interaction  of  a  mid-lati- 
tade  baroclinic  wave  with  several  different  mountain-valley 
configurations.  The  purpose  was  to  determine  the  contribu- 
tion of  a  well-developed  baroclinic  wave  to  the  initiation 
of  a  monsoon  surge.  It  was  found  that  these  waves  can  ini- 
tiate a  surge.  Such  surges  are  weak  and  limited  in  extent 
but  they  do  exhibit  the  characteristics  of  observed  monsoon 
surges.  The  latitudinal  location  of  the  east-west  mountain 
range  has  a  large  impact  on  the  strength  of  the  disturbance 
and  thus  on  the  initiation  of  the  surge.  However,  the  size 
of  the  valley  is  nDt  an  important  consideration.  The  surge 
is  strongest  at  the  lowest  levels  and  dissipates  rapidly 
with  height.  Overall  results  indicate  surge  forcing  by  the 
baroclinic  wave  is  possible  but  other  forcing  mechanisms  are 
necessary  in  crder  to  simulate  the  stronger,  more  extensive 
surges  observed   during  the  winter   monsoon. 
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I .        INTRODUCTION 

In  recent  years  meteorologists  have  given  considerable 
attention  to  the  Asian  winter  monsoon.  Described  by  Chang 
and  Lau  (1980)  as  "one  of  the  most  energetic  heat  engines  in 
driving  the  earth's  atmosphere",  the  winrer  monsoon  is  the 
mosx  distinctive  component  of  the  Northern  fiemispaere  winter 
circulation.  I-s  fundamental,  large  scale  structure  has 
long  been  observed,  however,  much  about  the  monsoon  remains 
unknown.  The  rcle  of  the  monsoon  in  the  ia-eraction  of  the 
mid-latitude  and  tropical  regions,  the  variability  of  the 
monsoon  and  rhe  forcing  cf  small  scale  monsoonal  variations 
are    all   areas   of   active   research. 

A  monsoon  may  be  defined  as  a  three  dimensional,  plane- 
tary scale  wind  regime  that  reverses  direction  seasonally. 
In  the  case  of  the  monsoon,  the  establishment  of  the  summer 
cyclone  and  winter  anticyclone  is  governed  by  land-sea  dif- 
ferences and  the  varying  position  of  the  sun.  Essentially, 
it  is  the  seasonal  replacement  of  one  persistant  circulation 
system  by  an  equally  persistant  but  opposite  system.  The 
development  of  these  circulations  closely  links  the  monsoon, 
particularily   that   of   east   Asia,    to   regional   topography. 
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The  long,  clear  nights  of  the  Sibarian  wintar  allow 
strong  radiaticnal  cooling  and  the  subsaquent  formation  of  a 
shallow,  cold  layar  of  surface  air.  Persistan-  cold  air 
advection  throughout  the  troposphere  further  contributes  to 
the  formation  of  a  surface  high.  The  winter  monsoon,  which 
dominates  the  entire  Asian  continent,  is  characterized  by 
this  intense,  cold  polar  anticyclone.  Tha  resulting  north- 
northeasterly  surface  flow  extends  over  all  of  eastern  Asia 
and  southward  across  the  equator.  Tha  Himalayan  Mountains 
restrict  the  southward  movement  of  tha  cold  lower  tropos- 
pheric  air  and  provide  an  effective  barrier  to  tha  mixing  of 
the  cold  continental  airmass  with  the  warm  tropical  airmass 
of  southern  Asia.  As  tha  anticyclone  becomes  established 
over  Siberia,  a  surface  cyclone  develops  over  tha  "maritime 
continent"  region  of  Indonesia  and  Malaysia  as  a  major  plan- 
etary scale  convective  area  is  shifted  eastward  from  its 
summer  position   near    India      (Ramage,    1971) . 

An  analysis  of  the  monthly  mean  velocity  potential  field 
at  200  mb  for  December  197a  by  Chang  and  Lau  (1980)  shows  a 
broad  region  of  divergence  over  the  maritime  continent. 
This  is  indicative  of  deep  rising  motions  and  large  scale 
overturning   of   the  atmosphere.         Their    analysis   also    shows   a 
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ragion  of  maximum  convergence  over  nortiiern  China  in  the 
vicinity  of  zhe  Siberian  anticyclone.  (Jlote  -hat  for  large 
scale  flow  the  divergence  pattern  is  approximarsl y  given  by 
the  velocity  potential.)  Upper-level  convergence  is  also 
evident  in  the  equatorial  west  Pacific  Dcean  and  over  the 
east  African  coast.  Although  much  wealcer,  these  latter  two 
areas  indicate  the  existence  of  a  i^allcer  circulation. 
Ascending  over  the  maritime  continent,  one  branch  flows 
westward  across  the  Indian  Ocean  and  Arabian  Sea,  while 
another  branch  flows  eastward  to  descend  over  the  western 
Pacific  Ocean. 

Thus  the  winter  tropospheric  circulation  in  this  region 
is  dominated  by  a  strong  local  Hadley  cell;  a  thermally 
direct  circulation  involving  descent  in  the  cold  air  over 
the  northern  anticyclone,  southward  surface  flow,  ascent  of 
warm  air  over  the  maritime  continent  and  nor-hward  return 
flow    at    200    mb. 

The  large-scale  components  of  the  monsoon  are  consistent 
and  identifiable  from  year  to  year,  but  they  are  not  static. 
The  primary  features  not  only  vary  in  strength  and  location 
on  an  annual  time  scale  but  there  are  small  time  scale  vari- 
ations in    the    velocity,       temperature      and    pressure    fields    as 
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well.  One  particular  event  that  can  be  linked  to  short  time 
scale  variations  is  the  monsoon  surge.  Also  known  as  "cold 
surges"  these  events  occur  a-  intervals  of  several  days  to 
two    weeks    during   the    winter   monsoon    (Lim   and   Chang,    1981). 

The  characteristics  and  affects  of  -he  monsoon  surge 
have  been  discussed  in  several  observational  srudies  at  both 
the  synoptic  and  planetary  scales  by  Danielson  and  Ho 
(1969),  Chang  et  al  (1979),  Murakami  (1979),  Chang  and  Lau 
(1980,1982)  and  Chang  et  al  (1983).  TtiB  following  sequence 
of    events    for    a   monsoon    surge   is    taken   from   these    studies. 

Prior  to  the  occurrence  of  a  monsoon  surge  there  is 
intense,  deep  cold  air  advection  over  northern  Asia  accompa- 
nied by  the  passage  of  a  mid-tropcspheric  trough  that  inten- 
sifies rapidly  over  northern  Japan.  Th=  resulting  increase 
in  baroclinicity  is  associated  with  ancicyclogenesis  over 
Siberia  and  cyclogenesis  in  the  strong  baroclinic  zone  off 
the  east  China  coast.  The  east  Asian  Hadley  circulation  is 
strengthened  through  increased  subsidence  as  the  anticyclone 
intensifies.  The  east  Asian  subtropical  jet  accelerates 
near  Japan  as  the  intensification  of  the  Hadley  circulation 
increases  upper-level  ageostrophic  flow.  The  west  Asian  jet 
streak  varies  inversely  with  the  east  Asian  jet  and  so 
decreases    in    strength   over  Afghanistan   and   Pakistan. 
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The  monsoon  surge  is  initiated  as  the  upper-level  trough 
moves  eastward  and  the  pressure  gradiaar  across  the  easr 
China  coast  -ighteas.  The  normal  geosxrophic  balance  of  rhe 
anticyclonic  flew  is  broken  since  the  pressure  gradient 
force  is  now  greater  than  the  Coriolis  force  and  a  rapid 
or oss-iso baric  flow  from  high  to  low  pressure  begins.  This 
accelerating  flow  towards  lower  pressures  at  lower  tropos- 
pheric  levels  is  the  monsoon  surge.  Following  rhe  surge,  an 
increase  in  the  northeasterly  flow  can  be  traced  to  the 
equatorial   regions  within    12-2U    h.  Within    one   to  two   days 

convection  in  the  sguatorial  region  intensifies  and  sustains 
the    already      enhanced  Hadley   circulation.  The   upper-level 

outflow  from  the  equatorial  South  China  Sea  accelerates  not 
only  northward,  but  within  one  day  of  the  surge  accelerates 
east  and  west  along  the  equator,  thus  strengthening  the 
Walker  circulation.  The  surge  ends  as  the  mid-latitude 
trough/ridge  pattern  moves  sufficiently  eastward  so  that 
subsidence   no    longer    occurs  over    the   anticyclone. 

Lim  and  Chang  (1981)  and  Chang,  Chen  and  Millard  (1983) 
suggest  that  the  monsoon  surge  is  characterized  by  two 
stages  separated  by  an  interval  of  a  few  hours  to  one  day, 
depending    on      the  location  of      the   observing      station.         The 
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first  stage  is  the  pressure  surge,  cr  L=ading  edge  of  the 
downward  accelerating  air,  with  a  propagation  speed  charac- 
teristic cf  gravity  waves.  The  second  stage,  which  loves  at 
advective  time  scales,  is  due  to  a  frontal  passage  and  is 
characterized   by   a   sharp    decrease   in   the   surface   dewpoint. 

Few  studies  have  been  done  to  determine  rhe  cause  or 
nature  of  the  forcing  that  precipitates  a  monsoon  surge. 
Lim  and  Chang  (1981)  used  the  linearized  shallow-water  equa- 
tions en  an  equatorial  beta-plane  in  an  analytical  study  of 
monsoon  response  to  a  transient  mass  source  that  simulated 
intense    anticyclogenesis.  Their   results    showed     that    such 

raid-latitude  forcing  excites  mainly  Rossby  waves  after  an 
initial  period  of  gravity  wave  type  motions  associated  with 
strong  northerly  winds.  They  suggested  that  pressure  forc- 
ing, using  equatorial  beta-'plane  dynamics,  can  explain  the 
large  features  of  synoptic  scale  motions  in  the  eastern  mon- 
soon region.  However,  they  do  not  consider  other  effects 
such  as  orography  or  a  variable  mean  zonal  wind.  Nor  does 
their  method  account  for  the  shallow  vertical  structure  of 
monsoon    surges   or  upper    tropospheric  response   to   the    surge. 

The  purpose  of  this  study  is  to  examine  the  interaction 
of    topography,    specifically  a   high    east-west    mountain   range. 
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with  a  mid-iatituds  baroclinic  wave.  Ths  purpose  will  be  to 
determine  the  contribution  of  the  baroclinic  wave  zo  the 
initiation  of  monsoon  surges.  A  global,  six-level  primitive 
equation  model  will  be  used  to  simulare  rhe  meteorological 
conditions  prior  and  subsequent  to  surgs  initiation.  Ini- 
tial conditions  will  include  a  strong  baroclinic  jet  north 
of  the  mountain  range  and  a  small  disturbance  in  the  mean 
current;  the  disturbance  will  grow  through  baroclinic  insta- 
bility. Several  mountain-valley  configurations  will  be  used 
to  study  the  surge  and  numerical  solutions  will  be  compared 
with    observed   monsoonal    surges. 
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II.        DISCRIPTIDN    OF    MODEL 

The  numerical  aodel  used  in  these  experiments  is  a  ver- 
sion of  the  UCLA  potential  enstrophy  conserving,  general 
circulation  model  for  an  adiabatic,  dry  and  frictionless 
atmosphere.  The  governing  equations  for  this  model  are  the 
primitive      equations:  the   horizontal     momentum      equation, 

thermodynamic  energy  equation,  surface  pressure  tendency 
(continuity)  equation,  and  the  hydrostatic  equation.  These 
are    given    in    Appendix   A.  The    design   and    dynamic   processes 

of  the  model  are  fully  discussed  by  Araiciwa  and  Lamb  (1977, 
1981).  A  summary  of  the  version  used  for  this  experiment  is 
presented    here. 

A.        THE    HCEIZCNTAL   GRID 

The  prognostic  variables,  wind(U,7)  ,  temperature  (T) , 
and  terrain  pressure (n),  are  staggered  horizontally  accord- 
ing to  Arakawa's  Scheme  'C  (Fig.  1).  This  grid  irrange- 
ment,  although  using  more  computer  time  than  an  unstaggered 
grid,  gives  a  better  structure  for  the  shorter  waves  that 
arise  due  to  geostrophic  adjustment,  nDcilinear  interaction 
or    mountain   effects.  Resolution    is   2.5°    long      by   'i.O^'    lat 
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using  a  36  x  46  spharical  grid.  Th*  coapatational  domain  is 
a  9C  sector  of  the  globe  wi-h  cyclic  contmaity  inipcssd  at 
tha    east    and   west  boundaries. 

In  spherical  coordinates  tha  poles  are  singular  points. 
Since  velocity  components  cannot  be  defined  at  the  poles, 
they  are  taken  as  n-points.  Due  to  the  convergence  of  tha 
meridians  towards  the  poles  a  very  short  time  step  would 
normally  be  required  to  insure  computatioaai  stability.  To 
allow  a  longer  time  step,  longitudinal  averaging  is  done  for 
selected  terms  in  the  equations  that  govern  gravity  wave 
propagation.  Tc  do  this  the  zonal  pressure  gradient  and 
zonal  mass  flux  are  expanded  into  a  Fourier  series  around 
each  latitude  circle,  then  the  amplitude  of  each  wave  compo- 
nent is  reduced  by  the  factor  S.  This  technique  smooths  the 
gradients  of  the  meteorological  fields  rather  than  their 
Fourier   expansions.  The  smoothing      operator    S      is    applied 

only    at    higher   latitudes    (where    S<1)    and   is   given    by: 

q   =  A^  c;os(|) 


AcJ)   sinCnAA/2) 


(Symbol    definitions    are    given   in    Appendix    3| 
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B.        THE    VEETICAL    GRID 

The    equations  used      in  this    modsl    are      writran    asing   the 
vertical    sigma   coordinate: 


p  -  p 
a  =              t  where     n  =   P     -  P^ 

— — s  t 


Ps    is     the   surface  pressure     and    P^      is   the   pressure      at    the 
tropopause . 

The  vertical  domain  extends  frora  the  earth's  surface  to 
200  mb.  Variables  are  staggered  vertically  (Fig.  2)  in  six 
layers  equally  spaced  in  sigma.  The  lower  boundary,  ct  =  1, 
follows  topography,  where  the  surface  pressure  is  defined. 
At  both  the  upper  and  lower  boundaries  vertical  velocity,  cr, 
vanishes.  The  prognostic  variables  (u,v,r,)  are  computed  at 
mid-layer,  while  geopotential ,  |,  vertical  velocity,  a,  and 
terrain  pressure,  n  are  computed  diagnostically  on  sigma 
le  vels. 

C.        FINITE    DIFFERENCING    SCHEMES 

The  time  differencing  methods  used  in  this  model  are  the 
fourth-order  space  differencing  scheme  and  the  !latsuno 
scheme    (also    called    the    Euler-backward   scheme) .      One    Matsuno 
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time   step    is   performed   followed  by    four   centered   time    s-eps. 

The   Matsuno      scheme    is   used      to  elimina-e      -he   computa- ional 

mode    of    the   centered    scheme  and  because   it    selectively    damps 

high    frequency   waves.      The  time  step   used    here   is    360   s. 
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Figure   1.       Horizontal   Distribution   of   Variables    on   a 
Spherical    Grid. 
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'igure  2.      Distribution    of  Prognostic   Variables    in   -he 

Vertical.      Solid  lines  are   a-levels,    dashed   lines 
are   at  aid-layer. 
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III.        XODEL    INITIALIZAIION 

The  model  uses  an  analytical  initialization  schsmo 
developed  by  Major  J.  L.  Ha  yes, i  OSAF.  A  small  perturba- 
tion is  superimposed  on  a  horizontally  and  vertically  vary- 
ing mean  current,  then  allowed  to  grow  through  baroclinic 
instability.  The  scheme  was  derived  specifically  for  flow 
over  a  flat  earth;  initialization  over  high  terrain  requires 
an    adjustment    to   the    scheme  and    is    discussed    separately. 

A.        INITIALIZATION   OF    THE    MEAN    FLOW 

A  baroclinic  Bickley  jet  is  situated  at  450N,  where  it 
reaches  a  maximum  velocity  of  46  m/s  at  200  mb.  The  initial 
mean  surface  flow  has  the  same  latitudinal  variation  as  the 
jet;  it  reaches  a  maximum  of  5  m/s  at  ^S'^N  and  decreases  to 
zero  at  the  poles.  The  v-coraponent  of  the  wind  is  initially 
set  to  zero.  The  surface  wind  is  described  by  the  Bickley 
jet    profile: 

U    C(j))      =     u     sech2(YCc^-<i)o)  ) 

S  IB 

V    C<j>)       =0 

s 


1    Personal   communication,    present   affilitation    Detachment    1, 
HQ    AWS. 
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The    upper-level   wind    has    the  same   latitudinal   profils 


U^((i),P(a,({)))      =     u     sech2(YC(J)- 


V^(<J»,P(a,c|)))      =     0 


The  lasr  term  en  the  right  produces  the  var-ical  shear.  The 
initial  wind  at  any  latitude  and  levsl  is  calculated  by 
adding  the   upper-level  and  surface   winds. 

Since  the  mean  flow  is  in  geosrrophic  balance,  the  sur- 
face pressure  is  calculated  by  intagrating  the  geostrophic 
wind   equation: 


p   ((^)      =     Pg(<j>o) 


(■ 


®^P  /"        "    s   tanh(Y(4>-4> 


RTqY 


u)A 


The    value      of    T, ,         288K,      is      taken    from    the      NACA   standard 
atmosphere   using    a  standard  surface   pressure    of    1013.25    mb. 

Temperature,      except      at   the    surface,       is      determined    by 
integrating   the   geostrophic  thermal    wind   aquation: 


T(p,<l))       =      T(p,4)o)    - 


U      a  2Ji 
m 


R  InCP/P 


t1 

m    J 


<1>0 


Integrals    are   computed   using   the    Simpson   Rule   approximation. 
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The  vertically  varying  term  T(p,0.)   is  obtained  from  lae 
NACA  standard  atmosphere: 

TCp,^0)   =   T(po,*o)   -   ^Z 
where:   Z  = 


■  ?  ('  -  (fen 


The   constants   are  defined   in    Appendix    B. 

B.        INITIALIZATICN   OVER    HIGH    TERRAIN 

The  mountains  in  this  model  are  intended  to  be  a  very 
simple  representation  of  the  Himalayan  Range.  The  mountains 
extend  east  tc  west  across  the  grid  section  wizh  a  maximum 
elevation      of    3000      m.  The   longitudinal      span      is    Ib'^      or 

approximately    1600   km. 

The  mountains  were  initialized  at  full  height  without 
any  growth  period.  The  equations  used  to  produce  the  moan-- 
tains  are    relatively    simple  and    produce    i    Scnife-edged    range: 


.25    ^  at    1  (^'<^i\    =   2A<^ 

$      =      Ci((!))$        =  I    .75    ^^        at    \(p"^i\    =   2A4) 

$  at   (})   =   !^1 


The  immediate   initialization  of   the  terrain   does  lengthen 
the  gecstrophic  adjustment  time.    However,  the  model  is  run 
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for  72  h  prior  to  creating  the  valley  aad  it  is  reasonabla 
to  assume  the  longer  adjustment  period  will  have  no  serious 
effect  on  the  attempted  simulation  of  the  monsoon  surge. 
Using  "he  mean  flow  initialization  scheme  without  a  mountain 
growth  period  resulted  in  a  large  scale,  north-south  oscil- 
lation in  the  pressure  field.  There  was  little  evidence  of 
damping  during  a  95-h  model  run.  To  solve  rhis  problem  the 
surface  pressures  at  mountain  grid  points  were  adjus-ed  for 
height  above  mean  sea-level.  Jsing  the  hydrostatic  equation 
and  the  adiatatic  lapse  rate  for  the  NACA  standard  atmos- 
phere, the  following  equation  was  developed  for  a  new  (moun- 
tain) surface  pressure  in  terms  of  the  old  (flat  earth) 
surface    pressure: 


P    (new)      =      P    ( 
s  s 


ox.  (.  -  I^) 


g/RT 


The  new  surface  pressure  is  used  to  re-calculate  wind  and 
temperature  values  at  each  mountain  grid  point.  No  oscilla- 
tion was  evident  when  the  adjusted  values  were  used  in  the 
initialization. 

As  the    monsoon   surge    is   expected   to    occur    almost    immedi- 
ately  following    the    creation   of    the    valley,      the   geostrophic 
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adjustmant  time  cannot  be  neglactsd  at  this  point,  Thars- 
fors,  the  valley  is  "built"  over  a  12-h  period  zo  minimize 
the   time      required   for   geostrophic   adjustment.  The    method 

used  is  based  on  the  scheme  developed  by  Hayes  (previously 
cited)  for  raising  terrain  height  during  a  model  run.  To 
build  the    valley   however,    the   method   is   reversed: 


*     = 


ClC*)*^  fl   -    C2(T)G(X)1 


The    function   C,(t)    is    defined: 


sm' 


C2(T)       = 


m) 


1  f 


T^T.    and   T>T- 
1  f 


The    function    G(X)    is    defined: 


G(X)   = 


X-X1-2AX 
-2AA 

A2+2AX-X 
2AX 

0 


Xl<X<X2 
(X-2AX)iX4Xi 

X2^X^a2+2AX) 

elsewhere 


X,  and  X^  are  the  west  and  east  edges  of  the  valley  respec- 
tively. Where  no  valley  is  built,  G  (X)  =1.  The  valley 
walls  are  sloped  over  two  grid-lengths  (2AX)  ,  and  the  valley 
floor  is    at   mean   sea- level    (0    m)  . 
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C.        TH2    PERTURBATION 

The  perturbation  is  a  standard  barotropic  disrarbancs. 
It  consists  cf  a  weak  longitudinal  sins  ^a  vs  with  rnaximum 
latitudinal  amplitude  at  450N.  The  vertical  aaplinude  is 
one-tenth  the  earth's  radius.  The  perturbation  is  described 
by: 

*•      =      fo(A  sin(nX)    sin^  (2c^)) 

The    following   perturbation   quantities      are    derived   using   the 
hydrostatic   equation    and    geostrophic   wind    equation: 

P»      =     _  ^'      =     -  ^"4  sin(nA)    sin^  (2(j>)) 
RT  RT  V  / 

IS*'  1 /  \ 

U'      =     -  ^  TT       =      -  -{A  sin(nA)    2sin(4(j))) 

afo  5(j)  a  V  / 


afo  cos      5  (|i  a  cosc^ 

T'       =      0 


1 [nA  cos(nX)    sin^(2(j))j 


The  total  fields  are  calculated  by  adding  the  perturbation 
quantities  to  the  aean  current.  The  total  fields  will  then 
be  integrated  in  time  to  arrive  at  the  required  forecast 
values. 
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I*^-       ANALYSIS    OF   lESULTS 

A.        EXPERIMENTS    1    S    2:       CONTROL    RUNS 

Two  central  runs  were  made  for  comparison  with  subse- 
quent model  runs.  The  first  run  uses  a  flat  earth  while  the 
second  includes  a  3000-m  east-west  mountain  range  cen-cered 
at  30°N.  Both  ruas  are  initialized  with  the  barcclinically 
unstable  mean  flow  and  the  barotropic  disturbance  discussed 
in  Chapter  3.  Model  times  given  are  based  on  the  cime  (h) 
at    which   the   disturbance    is   introduced   into   the    mean    flow. 

The    first   run  follows   the   development    of    the   disturbance 
independent  of   terrain   effects.  Sea-lsvel    pressure    fields 

for  60,  72,  84  and  96  h  are  shown  schemet ically  in  Figs.  3 
and  4.  Potential  temperature  fields  for  the  same  times  are 
shown  in  Figs.  5  and  6,  Since  a  pie-shaped  section  of  the 
globe  is  represented  in  cartesian  coordinates  there  is  some 
distortion  of  the  fields.  This  is  primarily  evident  in  the 
tight  north-south  pressure  gradients  wnirh,  for  comparable 
wind    velocities,    appear    weaker   in   the    east-west    direction. 

The   disturbance    is   strictly    limited   to    the    mid-latitudes 
with    a   very   flat    pressure    gradisnt   in    both    the    polar    and   the 
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tropical  regions.  The  disturbance  moves  eastward  st  a  con- 
stant velccity  of  approximately  8.5°  svery  twelve  hours. 
Wave  development  is  realistic  wi^h  a  995  mb  low  and  a  1023 
mb    high    developing  after    78  h. 

Development  of  the  disturbance  in  tha  wind  and  tempera- 
ture fields  closely  parallels  t ha-  of  the  pressure  field, 
although  the  temperature  and  potential  taaperature  data  dis- 
play far  more  short  wave  noise.  As  no  surface  frictional 
dissipation,  diffusion,  or  smoothing  -achniques  are  used  in 
the  model,  some  noise  in  the  resulrs  is  sxpec-ced.  The  large 
scale  temperature  features  remain  discernable  throughout  xhe 
six  day  (144  h)  run  but  the  the  small  scale  features  are 
masked  by  noise  after  126  h,  sspecially  ar  higher  levels. 
(Temperature  and  wind  velocities  are  plottad  on  sigma  sur- 
faces; • higher  levels*  refers  to  increased  height,  or  lower 
sigma  values. ) 

After  96  h  particularily  strong  cyclogenesis  begins  to 
occur.  At  114  h  the  cyclones  intensify  rapidly  and  the 
ridge  is  drastically  weakened.  While  diffusion  or  smoothing 
applied  to  the  fislds  would  decrease  the  strong  gradients 
caused  by  cyclogenesis,  they  would  also  wsaken  the  gradients 
of    pressure   and   temperature  which   could  contribute   to    possi- 
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bl9  mcnsoon  surges.  For  this  rsason  neither  mathod  is  used. 
A  linear  drag  term  was  applied  tc  the  i#inds  in  the  lowes- 
level,  however  it  had  no  discernable  affect  on  the  strong 
cyclogenesis. 

The  mountain  range  used  in  second  control  run  is  tha 
basis  for      all   subsequent   terrain   configurations.  Fig.      7 

gives  a  three-dimensional  reprasenta-cion  of  the  range. 
(Note  the  vertical  scale  is  greatly  exaggerated.)  The  data 
from  this  run  are  used  to  determine  the  time  and  location  at 
which  the    valley     is    to    be  created.  rha    reguiremants   are: 

(1)  a  high  pressure  center  to  the  northeast  (to  be  consis- 
tent with  the  Siberian  anticyclone)  and  (2)  a  sea-level 
pressure  difference  of  25  ab  between  zha  high  and  low  ores- 
sure  centers.  The  latter  requirement  ensures  an  adequate 
northwesterly   current   ro    initiate   a    monsoonai   surge. 

Sea-level  pressure  and  potential  temperat-ure  fields  for 
60,  72,  84  and  96  h  are  shown  in  Figs.  3  through  11.  East- 
ward movement  of  the  disturbance  is  again  S-S**  per  twelve 
hours.  The  mosi:  noticeable  difference  between  this  experi- 
ment and  the  first  is  a  series  of  weak  highs  and  lows  that 
have  developed  along  the  northern  slopa  of  the  mountains. 
These  disturbances  migrate  eastward  at  a  spaed  comparable  to 
the    mid-latitude    wave. 


31 


Development  of  the  wave  is  slower  than  in  th?  previous 
case.  At  78  h  a  999  mb  low  an d  a  1019  ib  high  have  devel- 
oped, giving  a  pressure  difference  of  20  ib  compared  to  28 
mb    earlier.  This    represents  an      overall    weakening      of   the 

disturbance  and  is  not  due  to  the  weakening  of  either  the 
high   or   low  centers    individually.  It   appears   th^   aountain 

range  may  have  a  stabilizing  effect  on  the  growing  distur- 
bance. The  required  25  mb  pressure  difference  is  reached  at 
84  h.  Since  a  twelve  hour  period  is  required  to  create  the 
valley,  growth  will  be  initiated  at  72  h.  model  time  in  sub- 
sequent runs.  Siace  the  valley  must  be  located  east  of  the 
anticyclone  in  order  to  take  advantage  of  the  associated 
northwest  flow,  the  western  edge  of  ta=  valley  is  placed 
directly  south  of  the  leading  portion  of  northwesterly  winds 
at    Sa   h. 

Development  of  the  disturbance  in  th3  potential  temperar 
ture  and  wind  fields  is  very  similar  to  the  first  run. 
There  is  a  maximum  in  potential  temperatura  over  the  moun- 
tains but    this    is  simply    due   to    the   increased  elevation. 

The  small  scale  noise  apparent  in  th=  potential  tempera- 
ture fields  of  the  first  run  also  appear  in  this  experiment 
although    at    a      later    time,      possibly   dua    to      the    oiuch    slower 
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dsvelopment  of  the  cyclone-anticyclone  system.  The  strong 
cyclogenesis  and  intense  pressure  gradients  also  appear 
later  in  the  second  model  run.  Wave  devsiopment  is  approxi- 
mately six  hours  slower  when  the  mounraias  are  included  and 
this  corresponds  to  the  time  lag  observed  in  the  cyclogene- 
sis   and  noise   growth. 

It  should  be  no'ed  that  potential  temperature  is  used  in 
the  analysis  of  model  results  since  irs  fields  will  no':: 
include  effects  of  the  adiabatic  processes  occurring  in  flow 
over    mountains. 

B,        EXPERIMENT    3 

The  initialization  for  this  experiment,  and  all  subse- 
quent experiments,  is  provided  by  the  sama  analytic  scheme 
used  in  -he  control  runs.  Terrain  consists  of  an  east-wes* 
mountain  range  centered  at  SO^'N.  The  valley,  built  between 
72  and  84  h,  is  twelve  gridlengths  wid=  or  nearly  one-third 
the  length  of  the  mountain  range  (Fig.  12)  .  The  critical 
time  period  for  the  occurrence  of  a  surge  will  be  the  first 
24  h  after  the  valley  is  completed  (84-108  h)  .  It  is  during 
this  period  that  the  anticyclone  is  located  northwest  of  the 
valley. 


33 


Development  of  the  disturbance  through  72  h  is  described 
by  ccntrol  run  #2.  The  sea-level  pressure  field  3.1  84  h 
shows  a  trough  extending  across  the  valley  and  along  the 
northern  slope  of  the  eastern  range.  kt  the  sane  time  a 
small,  high  pressure  ridge  can  be  seen  in  the  northwest 
region  of  the  valley.  For  the  next  twelve  hours  -he  ridge 
continues  to  build  southward  until,  at  96  h,  i-  extends  the 
length  of  the  western  valley  wall  splitting  -he  trough  into 
two    separate      low  centers.  The   pressure      gradient    in      the 

tropical  region  south  of  the  valley  has  also  increased  (Fig. 
13) .  As  the  ridge  continues  to  expand  southward,  the  east- 
ern portion  of  the  trough  develops  a  very  strong,  tight 
pressure  gradient  in  the  valley  and  along  the  nor-hern  moun- 
tain slopes  (Fig.  14),  The  ridge  continues  to  s-rengthen 
and  remains  located  along  the  western  wall  through  114  h. 
Af-er  this  time  the  anticyclone  moves  eastward  to  a  position 
due    north    of   the    valley    and  the    ridge    begins    to    weaken. 

Also  noted  in  the  results  is  an  oscillation  in  the 
intensity  of  the  pressure  field.  It:  begins  when  valley 
growth  is  initiated  and  has  a  period  of  i pprcximately  24  h. 
The  observed  oscillation  is  possibly  -he  result  of  gravity 
waves  set  up  by  the  relatively  sudden  (12  h)  creation  of  the 
valley. 
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The  second  stage  of  the  monsoon  surge  was  earlier 
described  as  being  characterized  by  a  sharp  drop  in  surface 
dewpoint.  Since  this  is  a  dry  model,  and  as  observational 
studies  (Chang  et  al,  1979,  Chu,  1978)  show  zhis  stage  can 
also  be  associated  with  a  drop  in  surfacs  temperatures,  the 
potential  temperature  fields  are  analyzed  for  possible  surge 
characteristics.  (Beyond  this  point  •  t2apera"cure'  refers  to 
potential    temperature   unless   otherwise   sta-ed.) 

In  agreement  with  the  development  of  the  surface  pres- 
sure ridge,  a  weak  tongue  of  cold  air  has  formed  along  the 
western  valley  wall  at  84  h.  However,  no  similar  gradient 
exists  along    the    eastern    wall    (Fig.    15). 

Twelve  hours  after  the  valley  is  completed,  the  cold  air 
has  reached  the  southern  valley  boundary.  Minimum  tempera- 
ture at  240N  has  dropped  from  31  OK  at  84  h  to  303K  at  96  h. 
A  tight  temperature  gradient  exists  along  the  western  side 
of  the  cold  air  intrusion  and  also  along  the  northern  moun- 
tain slopes  where  temperatures  have  remaiaed  relatively  con- 
stant. The  cold  air  tongue  can  easily  by  seen  at  the  two 
lowest  levels  (Fig.  16)  .  At  the  mid-lavels  a  very  slight 
tightening  of  the  gradient  can  be  seen  over  the  western  val- 
ley   wall    when      compared    to  the   eastern    wall,         but    only    over 
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th9  northern  half  of  the  valley.  At  tha  highest  two  levels, 
gradients  over  the  -wo  valley  walls  are  nearly  identical  so 
that  the  apparent  changes  in  temperature  are  probably  due  to 
the  effect  of  terrain  variation  on  ths  slope  of  zhe  sigma 
surface. 

The  cold  air  continues  to  penetrate  southward  for  the 
next  twelve  hours  reaching  its  maximum  southward  extent  at 
108  h  (Fig.  17).  At  240N  the  minimum  usmperature  is  298K,  a 
-12K  change  frcm  84  h.  The  cold  air  does  not  penetrate  much 
beyond  the  southern  valley  border.  At  23<^N  -he  temperature 
change  is  approximately  -2K  in  24  h.  i^hils  southward  move- 
ment of  the  cold  air  ceases  at  108  h,  it  does  expand  east- 
ward, apparently  in  response  to  the  eastward  movement  of  the 
anticyclone      located    north     of  the      valley.  The   cold      air 

tongue  does  not  dissipate  as  rapidly  as  the  pressure  ridge, 
and  a  pool  of  relatively  cool  air  remains  in  the  western 
valley  region  until  the  arrival  of  the  next  mid-latitude 
cyclone. 

At  the  higher  levels,  sigma  =  7/12  and  9/12,  a  small 
tongue  of  cold  air  forms  between  78  and  34  h  above  the  west- 
ern valley  wall.  The  southward  extent  is  slightly  less  than 
that    at   the   lowest  level.         As   the   low-level    cold   air    tongue 
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moves  rapidly  southward  aftar  8U  h,  tha  cold  air  at  zhese 
lavels  moves  very  slowly,  only  four  degraes  between  3a  and 
102  h  at  sigma  =  7/12.  Little,  if  any,  eastward  expansion 
is  seen.  Beyond  108  h  the  very  slight  gradients  are  diffi- 
cult to  discern,  primarily  due  to  the  increasing  small  scale 
noise.  At  sigma  =  9/12  the  cold  air  aora  closely  follows 
the  behavior  of  the  cold  air  tongue  at  tha  lowest  level 
showing  continued  south  and  eastward  expansion  on  the  same 
time   scale,    although    of    less   intensity. 

To  verify  that  the  changes  in  the  poten-ial  temperature 
fields  at  the  upper-levels  are  essentially  due  to  the  use  of 
the  sigma  coordinate  system,  the  temperarura  (T)  fialds  were 
analyzed  at  850,  700  and  500  mb  between  34  and  102  h.  At 
850  mb  the  intrusion  cf  ccld  air  can  easily  be  traced.  At 
700  mb  the  cold  air  tongue  is  not  discarnabla  although  this 
may  partially  be  due  -o  the  adiabatic  warming  effects  appar- 
ent at  both  the  700  and  850  mb  levels.  At  500  mb  some  small 
scale  noise  is  introduced  by  flow  over  and  around  tha  valley 
but  there  is  no  evidence  of  the  southward  movemen-  of  cold 
air.  Therefore  it  is  reasonable  to  assuoia  that  -he  changes 
in  the  potential  temperature  fields  at  -he  higher  levels  are 
due   to  the   effect  terrain    variation      has   on    the   slope    of   the 
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sigma  surface.  The   apparent   movement    of      colder   potential 

temperatures  southward,  then,  is  merely  the  result  of  the 
sigma  surface  being  at  a  lower  elevation  in  the  valley 
region.  To  scne  extent  this  effect  is  seen  near  the  eastern 
valley  wall   at  the  lower    levels    as    well. 

At  78  h,  when  the  valley  is  only  partially  completed,  a 
southward  deflection  of  the  northwesterly  winds  ahead  of  the 
ridge  is  readily  apparent.  By  84  h  the  valley  is  complete 
and  the  north  winds  are  rapidly  moving  southward  (Fig.  18). 
Along  the  valley  wall  the  leading  edge  of  the  wind  surge  is 
highly  ageostrophic.  A  strong  southwesterly  flow  has  set  up 
in  the  eastern  region  of  the  valley  presumaoly  in  response 
to    the      tight    pressure  gradient      that   exists    there.  It    is 

interesting  tc  note  the  small  cyclone  that  has  developed 
downwind  of  the  western  valley  wall.  It  appears  to  be 
slightly  better  defined  than  the  pre-existing  low  to  the 
northeast  and  may  result  from  the  combined  effects  of  tae 
shear  line  and  of  cyclonic  flow  induced  in  the  descending 
air  by  the  conservation  of  potential  vorticity.  The  cyclone 
dissipates  as  the  north  wind  continues  to  push  south  and 
eastward  into  the  valley.  The  shear  line  remains  well 
defined      between   the      opposing      northerly    and      southwesterly 


38 


flows.  Th€  strong,  cross-isobaric  soathwest  winds  in  the 
eastern  valley  coincide  with  a  very  tight  pressure  gradient 
and    may   represent  a    sloping  front   in  the   valley. 

At  96  h  the  north  wind  regime  extends  -he  length  of  the 
western  valley  wall  and  is  spreading  eastward  across  the 
valley    (Pig.       19)-  At    102    h  the   north      winds   have   reached 

their  maximum  extent;  the  wind  regime  occupies  fully  half 
the  valley  but  does  not  extend  south  of  it.  The  flow  along 
the  valley  wall  is  no  longer  ageostrophic.  The  north  winds 
entering  the  valley  between  84  and  102  a  are  iden-ifiable 
only  at  the  very  lowest  levels,  in  agreement  with  the  defi- 
nition of  a  monsoon  surge  as  a  low-level  phenomenon.  By  108 
h  the  anticyclone  has  moved  sufficiently  eastward  that  the 
light  northerly  winds  associated  with  the  ridge  line  have 
replaced  the  strong  northwest  flow  previously  located  north 
of  the  valley  (Fig.  20).  The  flow  into  the  valley  is 
greatly  reduced  and,  without  this  support,  the  surge  weakens 
and  dissipates.  The  eastern  region  continues  to  be  domi- 
nated  by    southwest   flow. 

It  is  clear  that  the  movement  of  the  north  winds  through 
the  valley  correlates  well  with  the  development  of  the  cold 
air    tongue   and   pressure    ridge   discussed   earlier. 


39 


In  summary,  several  observations  can  be  made  from  the 
analysis  of  these  fields.  First,  all  three  fields  shew  the 
very  rapid  southward  movement  of  a  cold,  high  pressure  air^- 
mass  along  the  western  valley  wall  between  84  and  96  h. 
Much  slower  southward  movement  continues  through  108  h. 
Second,  the  airmass  expands  eastward  across  the  valley  in 
response  to  the  continued  eastward  movsaent  of  -he  anticy- 
clone. And  third,  the  airmass  never  penetrates  the  tropical 
regions   but   remains    within  the   confines  of    -he   valley. 

C.        EXPERIMENT    4 

In   this   experiment   the  mountain      range    is   again  cen-ered 
at    30ON.  The   valley,       built   between      72    and    SU    h,         is    24 

grid-lengths  wide,  ie.  two-thirds  of  tae  mountain  range  in 
length.  The  terrain  then  is  actually  the  small  east-west 
oriented  mountain  range  displayed  in  Fig.  21.  As  in  Experi- 
ment #3,  the  critical  period  for  the  surge  is  between  82  and 
108  h  while  the  anticyclone  is  located  northwest  of  the 
valley. 

The    first    72      h    of   development    are      described    by   Control 
Run    #2.         As   in   the    previous     experiment   the    creation    of   the 
valley  causes      an  oscillation      in   the      pressure    (ie.         mass) 
field.        The   oscillation      in   this  experiment    is      larger    than 
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that  of  the  previous  case  since  the  reduction  of  terrain  is 
greater.  The  rapid  increase  in  cyclogsnasis  also  occurs 
earlier,  supporting  the  assumption  that  the  mountains  exert 
a   stabilizing   influence    on  the   disturban::s. 

At  8U  h  a  small  southward  extension  of  the  high  pressure 
ridge  can  be  seen  over  the  eastern  mountain  slopa  (western 
valley  wall).  A  low  is  centered  over  the  northern  valley 
with  a  weak  trough  extending  southwest  over  the  valley  wall. 
A  tight  gradient  develops  south  of  the  low,  extending  over 
the  eastern  valley  wall  and  continuing  along  the  northern 
mountain    slopes.  The    trough   is   strongest      at    78   h      and   is 

subseguently  pushed  southward  by  the  expanding  high  pressure 
ridge.  The  ridge  advances  southward  during  the  next  twelve 
hours  reaching  the  southern  valley  edge  at  96  h  (Fig.  22). 
Thus  far  movement  of  the  ridge  is  quita  similar  to  Experi- 
ment #3,  although  it  does  appear  to  be  slightly  weaker  and 
not    as  readily    discernable. 

Between  96  and  108  h  the  mid-latituda  ridge  weakens  con- 
siderably. Concurrently,  a  strong  pressura  gradient  builds 
along  the  western  edge  of  the  ridge,  in  part  due  to  a  trough 
extending  southward  from  the  intense  cyolone  located  to  the 
northwest     (Fig.    23).         A    similar    gradient    is    seen    during    the 
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same  period  in  Experiment  #3  but  is  not  nearly  as  intense. 
At  both  102  and  108  h  there  is  evidence  zha.z  -he  ridga 
extends    slightly   beyond    the  edge    of    the   valley. 

The  development  and  movement  of  the  ridge  agrees  very 
closely  with  that  of  the  earlier  experiment  using  a  smaller 
valley.  In  both  cases  the  ridge  begins  development  while 
the  valley  is  growing,  moves  southward  along  the  western 
wall  at  the  same  rate,  reaches  the  southern  edge  at  nhe  96  h 
point  and  weakens  when  the  anticyclone  moves  to  a  location 
due  north  of  the  valley.  There  is,  however,  a  slight  dif- 
ference in  intensity.  The  third  experiment  shows  a  well 
defined,  easily  discernable  ridge  from  34  to  108  h.  In  the 
present  experiment  the  gradient  along  the  western  side  of 
the  ridge  is  greater,  but  the  ridge  itself  is  not  as  well 
defined. 

A  weak  tongue  of  cold  air  can  be  sesn  along  the  western 
valley  wall  at  84  h  (Fig.  24)  .  An  increase  m  the  tempera- 
ture gradient  is  apparent  in  the  extreme  western  valley  at 
the  lower  three  levels  although  it  is  very  weak  at  the 
uppermost  of  these  (sigma  =  7/12).  These  levels  correspond 
to  pressure  heights  lower  than  600  mbs  in  a  standard  atmos- 
phere.       Through   the      next  twelve   hours  the      cold    air    tongue 
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moves  southward  and  strengthens.  3y  96  h  the  cold  air  has 
reached  the  southern  edge  of  ths  valley  (Fig.  25).  The  cold 
air  tongue  is  much  better  defined  than  ths  pressure  ridge. 
It  is  nearly  identical  to  the  cold  air  tongue  seen  in  the 
western  valley  for  Experiment  #3.  Tha  minimum  temperature 
at  2U«N  has  fallen  from  3  10K  at  84  h  zd  303K  at  96  h.  The 
temperature   change  is   primarily    limited  to    the   western    third 

of    the  valley,      temperatures   in      the   eastsrn    -wo-t.hirds   have 

« 
remained      relatively    constant      following      completion    of      the 

valley. 

The  cold  air  continues  to  move  south  and  eastward 
through  the  next  twelve  hcurs.  Maximum  southward  extent  is 
reached  at  138  h.  At  240N  the  minimum  tamperaturs  is  298K 
and  represents  the  same  change  in  temperature,  -12K,  as  was 
seen  in  Experiment  #3  (Fig.  26)  .  Mhils  southward  movement 
of  the  cold  air  has  ceased,  it  does  continue  to  expand  east- 
ward well  past  this  time.  The  primary  axis  of  the  cold  air 
tongue  remains  along  the  western  vallay  wall  and  retains  its 
identity  through  120  h.  A  pool  of  cool  air  remains  in  the 
western  valley  for  another  twelve  hours.  At  sigma  =  9/12 
the  cold  air  is  discernable  well  past  108  h.  At  higher  lev- 
els the  small  scale  noise  has  more  effect  and  the  cold  air 
is    difficult   to   trace   past  this    time. 
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The  upper-lsvels  display  the  same  behavior  seen  in 
Experiment  #3.  A  slight  increase  in  eastward  expansion  at 
these  levels   is   seen,    possibly   due    to    the    much   wider    valley. 

As  with  the  pressure  and  temperature  fields  the  develop- 
ment of  the  north  wind  surge  in  this  experiment  closely  fol- 
lows that  of  Experiment  #3.  At  84  h  the  norxh  winds  have 
reached  mid-valley  and  continue  to  spread  rapidly  southward. 
Again  the  leading  edge  of  the  wind  surge  is  ageost rophic.  A 
strong  southwesterly  flow  has  set  up  in  the  valley  appar- 
ently due  to  the  tight  pressure  gradient  located  there  (Fig. 
27).  The   small     cyclones  seen      in   Experiment      #3   are      not 

apparent  here  although  cyclonic  flow  exists,  possibly  due  to 
the  convergence  of  the  southwest  flow  from  the  tropics  and 
the  opposing  flow  out  of  the  mid-latitude  anticylones.  This 
cyclonic    flow   weakens   considerably   after   34   h. 

At  96  h  the  north  wind  regime  extends  almost  the  entire 
length  of  the  valley  and  is  spreading  eastward.  The  south- 
ward extent  of  the  north  winds  along  the  valley  wall  is 
identical  to  that  of  Sxperieratn  #3  (Fig.  28).  Further  east 
the  perimeter  of  the  north  wind  regime  is  nearly  two  degrees 
farther  south  of  that  reached  in  the  earlier  experiment. 
The    mid-latitude     ridge    weakens    during     the   next      six    hours. 
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and  vith  its  ccntinued  eastward  movemsat,  results  in  a 
decrease  in  the  northwesterly  winds  supporting  ths  surge. 
Thus  only  a  very  little  southward  movement  takes  place 
between  96  and  102  h,  the  time  of  maximum  southward  extent 
of  the  north  winds,  '  The  ageostrophic  component  of  the  flow 
weakens  rapidly  as  southward  movement  along  the  valley  wall 
slows.  The  winds  become  much  weaker  but  do  not  decrease  in 
coverage  through  108  h  (Fig.  29).  A  slight  eastward  move- 
ment is  seen  as  the  ridge  moves  past  thi  mountains  and  north 
of  the  valley.  After  108  h  the  north  winds  diminish  and 
move  eastward  in  response  to  the  weakening  and  continued 
eastward  movement  of  the  ridge.  After  120  h  the  wind  surge 
is    no   longer    discernable. 

The  surge  is  identifiable  only  at  the  two  lowest  levels. 
At  Sigma  =  9/12  the  surge  is  weaker  and  less  extensive  than 
at  the  lowest  level.  It  dissipates  very  rapidly  after  102  h 
as    the  ridge   at    this    level  moves    north   of    the   valley. 

In  this  experiment  a  slight  increase  in  the  eastward 
extent  of  the  north  winds  is  noted  when  compared  to  Experi- 
ment #3.  Since  topography  is  the  only  parameter  changed,  it 
is  reasonable  to  assume  the  increase  in  eastward  extent  is  a 
result  of    the      expanded    valley.         It   should      be    noted    though 
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that  while  tha  valley  size  has  double!,  the  corresponding 
increase  in  longitudinal  coverage  by  the  north  winds  is 
extremely    small. 

Several  observations  can  be  made  frDin  the  results  of 
Experiment    #4.  First,      all      three   fields      show   the      rapid 

southward  movement  of  the  cold,  high  pressure  airmass  along 
the  western  valley  wall  between  84  and  96  h.  Subsequent 
southward  expansion  is  much  slDwer.  Second,  the  airmass 
expands  eastward  across  the  valley  in  response  to  -he  east- 
ward movement  of  the  anticyclone.  And  third,  there  is  a 
tendancy  towards  cooling  and  rising  pressures  just  south  of 
the  valley  although  the  north  winds  remain  within  the  valley 
confines. 

D.        EXPERIMENT    5 

In  this  experiment  the  mountain  range  is  centered  at 
SS^'N,  eight  degrees  farther  north  than  in  previous  experi- 
ments. The  valley  is  twelve  gridlengths  wiie,  as  in  Experi- 
ment #3.  Terrain  is  shewn  in  Fig.  33.  Cyclogenesis  is 
considerably  slower  than  that  seen  in  either  of  Experiments 
3  or  4.  This  is  due  to  the  placement  of  the  northern  moun- 
tain slopes  at  the  same  latitude  as  the  maximum  amplitude  of 
the    disturbance.      Srowth    cf  the    perturbation    is   prevented   by 
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effectively   blocking      che    north^south    heat    flux.  For    this 

reason  the  initia-ion  of  valley  growth  was  delayed  six  hours 
to  allow  additional  time  for  development  of  the  disturbance. 
The  growth  period  for  the  valley  is  now  benween  78  and  90  h. 
As  with  the  previous  experiments  the  creation  of  the  valley 
causes  an  oscillation  in  the  pressure  field,  the  approximate 
amplitude    of  that  in    Experiment    #3. 

At  90  h  a  strong  Idw  is  located  in  che  northeast  of  the 
valley  trough  extending  southwest  over  the  valley  wall. 
This  results  in  a  tight  pressure  graiianz  in  the  eastern 
valley  and  along  the  northern  mountain  slopes  similar  to  the 
gradients  found  in  previous  experiments.  The  mid-la-itude 
anticyclone  is  neither  as  strong  nor  as  wrll  developed  as  in 
Experiments  3  and  U,  and  although  there  does  appear  to  be  a 
slight  southward  extension  of  the  high  pressure  along  the 
western  valley  wall,  it  is  very  weak.  During  the  next 
twelve  hours  a  ridge  of  high  pressure  builds  along  the  val- 
ley wall,  however  it  is  difficult  to  determine  whether  the 
high  pressure  originates  in  the  anticyclone  to  the  northwest 
or  in  the  tropical  regions  south  of  the  valley  (Fig.  31). 
As  the  analysis  of  the  wind  fields  will  later  show,  the  high 
pressure    air    does  originate  in   the   mid-latitude   anticyclone. 
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At  102  h  the  anticyclone  is  ill-defiaai  and  rhe  mountain 
region   zo    its   south    populated   by    small   high   and   low   centers. 

Between  102  and  114  h  the  ridge  moves  east  into  the  va.!- 
ley  and  strengthens,  this  is  par ticularily  noticeable  at  114 
h       (Fig.         32).  It   appears      the      anticyclone   has     slipped 

beneath  the  cyclone  and  is  re-forming  in  the  northern  val- 
ley. During  this  time  a  weak  trough  :;an  be  seen  forming 
over   the    northwestern   valley   wall. 

At  120  h  the  anticyclone  is  located  in  the  northeast 
valley  region  with  a  weak  ridge  extending  southwest  througn 
the  valley  and  a  stronger,  better  defined  ridge  to  the 
northwest.  Low  pressures  are  found  over  the  western  valley 
wall  due  to  a  weak  trough  forming  ahead  of  the  mid-latitude 
cyclone. 

In  earlier  experiments  the  high  pressure  ridge  developed 
and  strengthened  along  the  western  valley  wall  only  to  dis- 
sipate as  the  anticyclone  moved  to  a  position  northeast  of 
the   ridge.  In   this    experiment      the   ridge    is      weaker,      not 

well-defined  and  moves  eastward  with  the  anticyclone.  Since 
the  northern  valley  is  located  at  the  latitude  of  maximum 
amplitude  for  the  disturbance,  creating  the  valley  allows 
the    anticyclone      to    reform     in   this   area      as    the      weak    ridge 
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moves  eastward   into    ths    valley.  The   sams   reasoning    can   be 

used  to  explain  the  strong  development  experienced  b-j;  the 
low  during  the  valley  growth  period;  it  was  located  along 
the  northern  mountain  slopes  at  the  saae  longitude  as  the 
valley. 

A  weak  tongue  of  cold  air  can  be  seen  along  the  western 
valley  wall  at  90  h  (Fig.  33).  This  gradient  is  found  in 
the   lowest      three  levels,      weakening      with    height.  At    the 

upper-levels  (sigma  =  5/12,  3/12,  1/12)  nearly  identical 
gradients  exist  at  both  ends  of  the  valley  and  are  due  to 
the  effects  of  terrain  on  the  slope  of  sigma  surfaces.  The 
cold  air  tongue  strengthens  and  moves  rapidly  southward  in 
the  following  six  hours.  By  102  h  the  cold  air  has  reached 
the  southern  edge  of  the  valley  (Fig.  34|  .  The  development 
and  movement  of  the  cold  air  tongue  is  much  better  defined 
than  that  cf  the  pressure  ridge.  It  is  nearly  identical  in 
development  and  intensity  to  the  cold  air  tongues  seen  in 
Sxperiments  3  and  4.  The  minimum  temperature  at  32°N  (a 
comparable  latitude  to  240N  relative  to  the  position  of  the 
mountains)  is  295K  compared  to  306K  at  90  h.  This  is  a 
temperature  change  of  -11K  compared  to  -12K  and  -10K  for 
Experiments   3      and  4      respectively.         Again      the   temperature 
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change  is  primarily  limited  to  the  wsstarn  third  of  -he 
valley. 

The  cold  air  continues  to  move  south  and  eastward  during 
the  next  twelve  hours,  reaching  its  maximum  southward  extent 
at  1  14  h  (Fig,  35),.  At  32°N  the  minimum  -cemperature  is 
293K,  a  temperature  change  of  -13K  compared  to  -12K  in 
Experiments  3  and  4  for  a  2U-h  period  following  complerion 
of  the  valley.  Eastward  expansion  of  -hs  cold  air  continues 
past  114  h.  In  the  previous  experimeics  the  axis  of  the 
cold  air  tongue  remained  along  the  western  valley  wall.  In 
rhis  experiment  the  cold  air  moves  eastward  as  well,  parri- 
cularily  in  the  northwest  region.  This  is  probably  the 
result  of  the  anticyclone  reforming  in  the  northern  valley 
and  the  subsequent  eastward  movement  of  the  pressure  ridge. 
The  cold  air  tongue  does  not  move  as  far  eastward  as  the 
ridge  and  begins  to  weaken  when  the  anticyclone  moves  north- 
east  of    the   valley. 

The  upper-levels  display  the  same  behavior  seen  in 
Experiments  3  and  4.  A  very  slight  movement  of  the  northern 
por-icn  of  the  cold  air  axis  is  seen  at  sigma  =  9/12,  but  no 
similar    movement     is    seen    at      sigma    =   7/12.  The    increased 

small     scale  noise      at   these     levels    makes      analysis    of      the 
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temperature  fields  difficult  past  114  h.  The  highsst  levels 
(sigma  =  5/12  to  1/12)  show  no  sign  of  sDuthward  laovsment  of 
cold   air    at   any   time. 

As  in  the  previous  experiments  southward  deflection  of 
the  northwest  winds  into  the  valley  occurs  before  the  valley 
is  completed.  By  90  h  the  norrh  winds  have  reached  mid-val- 
ley with  a  shear  Line  developing  over  the  valley  floor  in 
advance  of  the  northerly  flow.  North  winds  along  the  valley 
wall  are  strongly  ageost rophic.  Again  a  strong  southwest 
flow  has  set  up  in  response  to  the  tight  pressure  gradient 
in  the  eastern  valley.  Although  the  valley  is  the  same  as 
that  used  in  Experiment  #3,  the  subsequent  pattern  of  devel- 
opment is  slightly  different.  At  96  h  the  northwest  winds 
reach  the  eastern  valley  wall,  whereas,  in  Experiment  *3  the 
strong  southwest  flow  in  this  region  prevents  the  northwest 
winds  from  spanning  the  valley.  The  wind  surge  is  actually 
strongest  at  this  time,  and  has  its  greatest  southward 
extent  (Fig.  36)  .  By  102  h  the  ridge  is  north  of  the  valley 
wall  and  the  winds  are  beginning  to  wealcen,  the  flow  is  no 
longer  ageostrophic.  The  north  winds  never  reach  the  south- 
ern edge  of  the  valley  (Fig.  3  7).  The  surge  continues  to 
weaken  rapidly      and    at      108  h      the   anticyclone      is   reforming 
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within  the  valley.  North  winds  exist  in  the  central  vallay 
region  but  dc  not  coincide  with  the  positon  of  the  cold, 
high  pressure  ridge.  Six  hours  later  the  anticyclone  is 
esxablished  over  the  north  central  vallay  and  no  evidence  of 
the    surge    exists    (Fig.      38).  It   is   clsar   that   positioning 

the  mcun-ain  range  at  380N  has  a  detriaanzal  effect  on  the 
development   of  the  surge. 

The  surge  is  identifiable  at  only  tha  "iwo  lowest  levels. 
At  Sigma  =  9/12  tha  surge  is  less  extensive  and  weaker  than 
at  the  lowest  level.  At  sig  ma  =  7/12  and  5/12  a  wealc 
cyclonic  turning  of  the  winds  is  apparent  but  no  dsfinite 
southward  flow  can  be  seen.  The  cyclonic  flow  disappears  as 
the    anticyclone   reforms    over   the    valley. 

In  summary  the  following  observations  can  be  made  from 
this  experiment.  First,  all  three  fields  show  the  rapid 
southward  movement  of  cold,  high  pressure  air  along  the 
western  valley  wall  during  tha  first  six  hours  following 
valley  completion.  Second,  the  movement  and  development  of 
the  cold  air  tcngue  is  much  better  defined  than  that  of  the 
pressure  ridge  and  wind  surge,  both  of  -rfnich  weaken  consid- 
erably after  102  h.  Third,  the  pressure  ridge  and  wind 
surge   axis      move    eastward      off   the      valley    wall      due    to      the 


52 


development   of   the  anticyclone   within      -h9    valley,      only   the 
northern    :nost    portion   of    the   cold   air    tongue    moves   eastward. 
And    fourth,      although   there  is  a    slight  oooling    and  pressure 
rise   south  of   the  valley,    the   north   winds    remain   well    within 
the    valley. 

E.        EXPERIMENTS    6    AND    7 

The  valley  was  included  as  part  of  the  initial  terrain 
configuration  for  these  two  experiments.  The  mountains  are 
centered  at  30<>N  and  the  valley  is  eigit  gridlengths  wide, 
approximately  one-quarter  the  length  of  the  mountain  range. 
This  is  the  original  mountain-valley  configuration  used  with 
the  model,  the  larger  valley  was  used  in  subsequent  runs  to 
avoid  interference  between  the  strong  southwesterly  flow  in 
the  eastern  valley  and  developing  soutiiward  flow  in  the 
west.  The  southward  movement  of  a  cold,  high  pressure  air- 
mass  occurs  with  both  valley  sizes,  however  the  airmass  is 
more  easily  seen,  and  so  more  easily  analyzed,  with  the 
larger  valley. 

Only  the  mean  flow  is  initialized  in  Experiment  #6.  The 
purpose  of  this  run  is  to  confirm  the  assumption  that  the 
monsoon  surge,  if  it  occurs,  is  due  to  baroclinic  wave  and 
not   a     result   of   the    interaction      of   the   mean   flow      with   the 


i.    -iT-  — 


erram. 


As  expected  no  "surge"  occurs.  The  current  tha-  dees 
form  in  the  valley  is  southw ast eriy,  flowing  from  high  to 
low  pressures.  The  strong  southwesterly  flow  seen  in  Exper- 
iments 3  through  5  is  in  part  due  to  the  mean  pressure  gra- 
dient that  exists  between  the  northern  and  southern  valley 
borders.  By  30  h  a  small  cyclone  has  formed  over  the  west- 
ern valley  wall,  probably  due  to  the  conservation  of  poten- 
tial vorticity.  This  indicates  that  the  small  cyclone  seen 
here  and  in  other  experiments  is  a  result  of  the  mean  flow 
interacting   with  the    terrain. 

Experiment  #7  is  initialized  with  the  perturbation. 
During  the  first  portion  of  the  experiment,  while  the  dis- 
turbance is  weak,  southwesterly  flow  forms  in  the  valley  in 
response  to  the  existing  north-south  pressure  gradient. 
Again  at  30  h  a  small  cyclone  has  formed  over  the  western 
valley  wall.  During  the  next  twelve  hours  the  center  of  the 
cyclone  moves  ever  tha  valley  floor  althougn  the  north  winds 
remain  located  over  the  valley  wall.  Ihs  mid-latitude  dis- 
turbance continues  to  strengthen  and  by  66  h  the  northerly 
winds  are  being  enhanced  by  the  northwest  flow  ahead  of  the 
approaching   ridge.  The   northerly    flow    strengthens      as   the 

ridge   moves    eastward,    at    96  h    it    extends    slightly    beyond   the 


54 


valley  edge,  however  this  is  its  niaximam  southward  extent. 
Subsequently  the  mid-latitude  ridge  axis  moves  north  of  zhz 
valley  and  the  supporting  northwest  flow  diminishes.  The 
flow  along  the  valley  wall  is  never  ageostrophic  during  this 
time. 

The  cyclone  is  readily  apparent  in  the  sea-level  pres- 
sure fields.  However,  net  until  the  ridge  approaches  the 
valley  does  a  high  pressure  ridge  build  along  the  valley 
wall  (72  h)  .  The  ridge  weakens  as  the  anticyclone  moves  to 
a    position   north   of    the    valley. 

The  potential  temperature  field  shows  very  little  cool- 
ing along  the  western  valley  wall  compared  to  previous 
experiments.  Minimum  temperature  along  2'4^N  at  90  h  is  299K 
compared  to  303K  at  66  h,  a  decrease  of  only  4K  in  24  h.  A 
very  small  temperature  change  when  compared  to  Experiments  3 
through    5. 
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Figure   7.      Model   Terrain,    Experiment   2. 
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Figure  8.       Sea-Levsl    Pressure,    60    and   72    h.    Ex?.    2 
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Figure    12.      ."lodal  Terrain,   Sxperiment    3, 
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Figure   14.      Sea-Level   Pressure,     108    h,    2xp.    3 
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Figure  21.      Mcdel  rsrrain,    Experimen-    fx 
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Figure   22.      Sea-Level   Pressure,     3u    and    96    h,    Sxp.    4. 
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Figure    23.      Sea-Level    Pressure,     108    h,    Zxd.    U 
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Figurs  26.   Potential  Tamperature,  108  h,    Exp.  U.   Sigma  = 
9/12,  11/12. 
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Fiaure   30.      Model  Terrain,   Experiment   5 
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Figure   33.      Potential  Temperature,    90   H,    Exp.    5.      Sigma    = 
9/12,    11/12. 
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V.        CON CLa  SIGNS 

The  interactioa  of  a  mid-latitude  baroclinic  w=v3  wi-ih 
several  different  mountain-valley  configurations  was  exam- 
ined. The  purpose  was  to  determine  the  contribution  made  by 
a  well-developed  baroclinic  wave  to  tha  initiation  of  mon- 
soon surges.  It  was  found  that  these  waves  can  initiate  a 
monsoon  surge  although  such  surges  are  waaic  and  limited  in 
extent. 

The  mountains  in  each  experiment  ars  an  east-west  range 
located  south  of  the  mid-latitude  wave.  As  was  first  noted 
in  Experiment  #2  the  mountains  appear  zo  have  a  stabilizing 
effect  on  the  growth  of  the  wave.  The  affect  is  especially 
evident  when  the  mountains  are  placed  farthar  north,  nearar 
the    center   of   the  disturbance.  The   reduced   growth    rate   is 

probably  caused  by  the  blocking  of  tha  v-componant  of  the 
wind  by  the  mountain  barrier  to  the  south.  This  would  slow 
the    wave      growth   because      the   energy      conversion   is      porpor- 


tional  to  V»e'.  Walker  (1982)  also  found  this  effect  with  a 
disturbance  moving  along  the  north  side  of  a  mountain  range. 
Previously,  several  observations  vara  made  about  tha 
nature  of  monsoon  surges.  Ihesa  include:  (1)  the  occurrence 
of    rapid,      ageostrophic    flow    from   high   to    low    pressure,       (2) 
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strengthening  of  the  north-northeasterly  winds,  (3)  a  sharp 
drop  in  surface  teaperature  and  (4)  -he  shallow  s-ructure  of 
the  surge.  Tc  seme  extent,  each  of  these  cnaracteris::  ics  is 
simulated    by   one   or    more    of   the    experimear  s. 

Experiment  #6,  in  which  only  the  zonal  mean  flow  inter- 
acts with  the  terrain,  shews  no  evidencs  of  a  surge.  All 
other  experiments,  in  which  the  perturbation  is  added  to  rae 
mean  flow,  do  produce  weak  monsoon  surges.  This  indicates 
the  eastward-moving  anticyclone  associa.ei  with  the  mid-la- 
titude baroclinic  wave  is  the  source  of  the  surges  occurring 
in    this    study. 

In  Experiments  3  and  4  the  mountains  are  centered  at 
30°N,  the  western  edge  of  each  valley  is  located  at  the  same 
point  although  the  fourth  valley  is  twice  the  size  of  the 
third.  In  each  case  a  cold,  high  pressure  airmass  moves 
southward  along  the  western  valley  wall  during  the  twelve 
hours  following  coapletion  of  the  valley.  The  leading  edge 
of    the  surge   is   ageostrophic. 

Figs.  39  and  40  show  the  change  of  potential  temperature 
with  time  in  the  vail  ay  for  Experiments  3  and  4.  The  pro- 
files extend  from  14on  to  46°N  along  the  western  valley 
floor    (55°      long)  ,      the    valley      is    located    between      22<5N   and 
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33<^N.  In  each  case  little,  if  any,  change  in  temp-ratur^ 
occurs  at  la^N  although  at  18°N  the  temperature  decreases  7K 
between  96  and  108  h.  In  the  valley  rhe  aaximura  tempera-rare 
change  occurs  earlier,  between  B'l  and  96  h,  and  is  much 
larger,      -16K   and  -17K  at    26°N.  The   surge    is    well   defined 

within  the  valley  but  dissipates  rapidly  south  of  22°N, 
This  is  indicative  of  the  weak  character  of  the  monsoon 
surges  in  this  study  since  the  cold  air  in  observed  cases 
has  been  traced  well  into  the  tropical  regions  while  in 
these  experiments  it  cannot  be  found  even  eight  degrees 
south  of   the   valley. 

As  can  be  seea  the  results  of  the  two  experiments  are 
nearly  identical,  indicating  the  size  of  the  valley  is  not 
an    important   consideration. 

The  oscillation  in  the  pressure  (mass)  field  is  appar- 
ently a  gravity  wave,  the  result  of  an  inward  rush  of  air 
filling  the  "vacuum"  produced  by  the  creation  of  the  valley. 
It  has  roughly  a  21-h  period,  which  may  be  related  to  the 
twelve  hour  period  in  which  the  valley  is  built,  and  it 
slowly  dampens  during  the  subsequent  60  li  of  the  experiment. 
It  has  no  observable  effect  on  the  pressure  gradients  and 
was  not  evident  in  either  the  potential  temperature  or  the 
wind    fields.      The  oscillation   also    occurs    in    Experiment    #5. 
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In  Experiment  #5  the  rapid  southward  movement  of  the 
cold  airmass  occurs  only  during  the  first  six  hours  after 
valley   completion.  This  is      due    to      the    placement      of   the 

mountains  at  38°N  resulting  in  a  much  weaker  disturbance. 
Fig.  41  shows  the  change  in  potential  teiperature  with  time 
in  the  valley.  The  valley  is  located  between  30°N  and  U6°N, 
the    profile   is   along    the    western    valley   floor. 

The  decrease  in  temperature  is  evident  in  -he  valley 
between  90  and  102  h.  As  in  Experiments  3  and  4  the  cold 
air  can  be  traced  only  four  degrees  south  of  zhe  valley. 
Unlike  the  two  previous  experiments  warming  occurs  during 
the  second  twelve  hour  period,  due  to  the  re-development  of 
the    anticyclone   within  the   valley.  The    southerly   winds    in 

the  western  region  of  the  anticyclone  prevent  the  flow  of 
cold  air  into  the  valley  and  effectively  "shu-  off"  the 
surge. 

The  monsoon  surges  initiated  in  these  experiments  are 
evident  only  in  the  lowest  sigma  levels,  at  sigma  =  11/12, 
9/12  and  7/12.  In  a  standard  atmosphere,  over  fla-  terrain, 
these  correspond  to  pressure  heights  less  than  675  ib.  The 
surge  is  strongest  at  the  lowest  level  and  dissipates  with 
height,         becoming   very       weak      at   sigma      =      7/12.        This      is 
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observed  in  both  the  wind  and  potential  tamperature  fields 
in  all  experiments  and  agrees  closely  with  the  definition  of 
the    monsoon   surge  as    a   low-level,    or   surface,    phenomenon. 

These    results     suggest  that    the      mid-latitude    baroclinic 
wave   can    initiate  monsoon   surges.  The    surges    are  weak    and 

limited  -o  the  valley  but  otherwise  exhibit  t.he  characteris- 
tics of  observed  aonsoon  surges.  The  weakness  of  these 
surges  indicates  o-cher  forcing  mechanisms  are  necessary  in 
order  to  adequately  model  the  monsoon  surge.  These  may 
include  a  strong,  stationary  high  located  northwest  of  the 
valley  similar  to  the  anticyclone  associated  with  the  winter 
monsoon.  Other  analytical  studies  indicate  that  an  upper- 
tropospheric,  planetary  scale  wave  number  three  is  necessary 
for    the   occurrence     of   strong    monsoon   surges.  Such    a    wave 

positioned  with  the  ridge  to  the  west  of  the  surface  anticy- 
clone would  result  in  strong  northwesterly  flow  over  the 
high  pressure  center.  This  would  initiate  cold  air  advec- 
tion  throughout  the  troposphere  and  increase  subsidence  over 
the  surface  anticyclone  leading  to  stronger  surges.  Weak 
monsoon  surges  were  associated  with  near  zonal  flow  at 
higher  altitudes,  this  describes  the  upper-level  flow  used 
in    these    experiments.      The  inclusion   of   one   or    both   of   these 
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phenomena    may     result   in      stronger,      more      extensive    monsoon 
surges  than   these  observed  in   this   study. 
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Figure  39.  Potential  Tamperatare  vs  Larituie,  Exp.  3. 
Vertical  lines  reoresent  valley  boundaries 
84  h,  (b)  96  h  and  (c)  108  ^i. 
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Figure  liO. 


Potential   remperature    vs   Latitude,    3xd.    4. 
Vartical    lines  represent   valley   boundaries, 
84    h,    (b)    96    h  and    ( c)     108    h. 
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Figure  41,  Potential  Temperature  vs  Latitude,  Exp.  5. 
Vertical  lines  represent  valley  boundaries 
90  h,  (b)  102  h  and  (c)  14  h. 
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APPENDIX  A 


MODEL  PRIMITIVE  EQUATIONS 


Equation  of  State: 


pa  =  RT 


Hydrostatic   Equation 


11=  - 


da 


TTO 


Horizontal   Momentum   Equations: 


n/dv\  +  fk  nv  +  V  (n$)  -  vn  5_($ff)     =    o 


\dt/ 


6a 


Continuity   Equation: 


5n  +  7  • (nv)  +  6  (na)  =    o 


St 


&a 
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l3t   Law   of   Thermodynamics: 


6    "(ncplne)    +  V«  (nvCpln9)    +  5_   (naCplne)      =     0 
5  So 
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APPENDIX    B 


LIST   OF    SYMBOLS 


a 
A 

Cp 

f 
fo 

F 

g 

n 
P 

Po 

Pm 

R 

T 

To 

a 
am 

V 

z 

a 

y 
r 
e 
\ 

IT 

n 

(T 
& 

T 

T. 
i, 

T, 

r 

<f> 


radius   of  the  earth 

amplitude  of  the    perturbation 

specific  heat   at    constant    pressure 

ccriclis  parameter 

value   of  coriolis  parameter   ar   45N 

force 

gravitat.ional  acceleration 

wave    number 

pressure  (P   =    ?  +   P'  ) 

surface    pressure    at   45N ,    1013.25    mb 

maximum   model  pressure    height,    125    mb 

gas   constant   for    air 

temperature 

surface   remperature   at    45M,    288K 

zonal   velocity   component      (U  =   U    +    U') 

maximum    velocity,    surface:    5    m/3,    upper:    65    m/s 

meridianal    velocity   component      (V    =    V    +   V') 

terrain   height 

specific  volume 

halfwidth  of    jet 

dry   adiabatic  lapse   rate,    6.5   K/!cm 

potential  temperature 

longitude 

3.14159 

normalized    pressure 

vertical  coordinate 

vertical  velocity   in   sigma   ccoridinate   system 

model   time 

model  tima    at   which  terrain   growth    begins 

model   time    at   which    terrain    reaohes    final    heigh' 

latitude 
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latitude  Df    j9t 

latitude  at    which  mountain   range    is    centered 

geopoter.tial,    also  current   terrain   heigh- 

maximum   terrain    height,    3000    m 

angular   velocity   of  the    earth 
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